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Abstract 

We analyze the vibrational resonance in the Duffing oscillator system in the presence of (i) a 
gamma distributed time-delayed feedback and (ii) integrative time-delayed (uniformly distributed 
time delays over a finite interval) feedback. Particularly, applying a theoretical procedure we obtain 
an expression for the response amplitude Q at the low-frequency of the driving biharmonic force. 
For both double-well potential and single-well potential cases we are able to identify the regions in 
parameter space where either (i) two resonances, (ii) a single resonance or (hi) no resonance occur. 
Theoretically predicted values of Q and the values of a control parameter at which resonance occurs 
are in good agreement with our numerical simulation. The analysis shows a strong influence of 
both types of time-delayed feedback on vibrational resonance. 

Accepted for publication in J. Applied Nonlinear Dynamics (in 2015) 

PACS numbers: 05.45.Ac, 05.45.Pq, 05.45.Xt 

Keywords: Vibrational resonance, Duffing oscillator, Distributed Delay feedback, biharmonic force. 


*Electronic address: rajasekar@cnld.bdu.ac.in 
^Electronic address: miguel.sanjuan@urjc.es 


1 



I. INTRODUCTION 


Time-delay is thought to be ubiquitous in many realistic models in physics, engineering 
and biology. When the state of a system at time t depends on its state values at previous 
past time, then a time-delayed feedback term is introduced in the evolution equation of 
the system. We can think about two main types of time-delays: discrete time-delays and 
distributed time-delays [l|. In the present work we are concerned with distributed time- 
delays. 

The general form of a distributed delay feedback term (DFT) in a dynamical system is 
given by 


DFT = 


G{T)x{t — t) dr. 


( 1 ) 


where a; is a state variable of the system and G{t) is a distributed delay kernel with G(r) > 0 
and G{t) dr = 1. The choice G as the Dirac-delta distribution 5{t — a) gives 


DFT = 


5{t — a)x{t — r) dr = x{t — a), 


( 2 ) 


that is, the delay time is a constant. The G{t) being 

Guir) = 


1/a, r G [0, a] 
0, otherwise 


( 3 ) 


leads to DFT = (1/a) x{t — r) dr which under a change of variable t — t = t' becomes 


DFT = 


a 


x{t') dt'. 


( 4 ) 


’ t—a 


The kernel in this case is uniformly distributed and the delay feedback term given by Eq 


is known as an integrative time-delay term 


a 


The parameter a represents the width and 


amplitude of the uniform distribution G^. The mean and variance of the time-delay in the 
distribution Gu is a/2 and a^/12, respectively, a can be treated as the response time of a 
^stem. The integrative time deW was earlier considered in the ‘integrate-and-hre’ models 
2] and self-organized criticality j^. This time-delay gives rise to amplitude death 0,1^, 7| 
in coupled oscillators and multiple resonance curves with and without hysteresis 

Another form of distributed delays which is shown to have a strong influence on the 
dynamics of a system is the gamma distribution of delays given by 


GJt) = 


^a^e 


Tip) 


( 5 ) 
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where a, p > 0 and r(p) is the Euler gamma function. For p = 1 the distribution Gg becomes 
an exponential distribution. For the gamma distribution (r) = p/a while = p/a^. 

The effect of distributed delays has been reported in a multitude of scientihc problems. 


Among them, we refer on stability 


10 | and emergence of chaos [n| in neural networks, 


Hopf bifurcation and stabilization of a fixed point in a discrete logistic model 12| , amplitude 
death in coupled oscillators [l3|, [l^, global asymptotic behavior of a chemostat model sys¬ 


tem 


15| . existence of wavefront solutions in reaction diffusion system s |l6l .ll7l|. bifurcation 


analysis in an epidemic model 18|, transmission dynamics of malaria jl9|], dynamics of car- 


20|, hrst passage time statistics |2l| and stochastic process in a Langevin 


following model 
equation [ 22 1 . 

Studies of different nonlinear phenomena in time-delayed systems have received a great 
interest in recent years. And the main goal of this paper is to investigate the role of gamma 
distributed time-delayed and integrative time-delayed feedbacks on the vibrational resonance 
produced in a Duffing oscillator, taken as the reference model system. Vibrational resonance 
is a resonant dynamics induced at the low-frequency u of the input periodic signal by a 
relatively high-frequency of the input signal 23l425|| . The effect of a single constant and 


multiple constant time-delays on vibrational resonance has been previously reported 


26- 


29|. When the delay time is not constant or not priorly known, then it is more realistic to 


consider distributed delays, and this is what we do here. 

The paper is organized as follows. In sec. 2 we consider the Duffing oscillator driven by a 
biharmonic force with two frequencies uj and D with D S> cn and with a gamma distributed 
time-delayed feedback term. We denote 7 as the strength of the feedback term. A gamma 
distributed time-delayed feedback is characterized by parameters 7 , a and p (refer Eq. ([5])). 
Assuming that the actual motion of the system consists of a slow motion and a fast motion 
and applying a theoretical method, we construct an analytical expression for the response 
amplitude Q (ratio of the amplitude of the slow motion of the output and the amplitude of 
the low-frequency component of the input biharmonic driving force). We analyze the double¬ 
well potential and the single-well potential cases separately. We treat the amplitude g of the 
high-frequency force as a control parameter. From the obtained theoretical expression of 
Q, we obtain the expression for g at which resonance occurs. For the double-well potential 
case, we identify the regions in the ( 7 ,p) parameter space where either two resonances or 
just one resonance occur. We point out the mechanism of occurrence of the resonance in 
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terms of the resonant frequency. In the single-well case, we show that at most only one 
resonance can occur. We identify the conditions on 7 and p for a single resonance and no 
resonance. We show that all the theoretical results are in good agreement with our numerical 
simulations. In sec. 3 we analyse the influence of the integrative time-delayed feedback term 
in the Duffing oscillator. The integrative time-delayed feedback term is characterized by 
the time-delay parameter a (Eq. (jl])) and the strength 7 of the feedback. For this feedback 
type, we are also able to obtain a theoretical expression for Q. We report the effect of the 
parameters 7 and a on vibrational resonance in detail for both, the double-well potential and 
the single-well potential system separately. The hnal section contains concluding remarks. 

II. DUFFING OSCILLATOR WITH A GAMMA DISTRIBUTED TIME- 
DELAYED FEEDBACK 

The dynamics of the Duffing oscillator subjected to a biharmonic force and a gamma 
distributed time-delayed feedback term is governed by the equation of motion 

X + dx + UqX + (3x^ F{t, x{t — t)) = f cos ut + g cos fit, ( 6 ) 

where D S> a; and F is the gamma distributed time-delayed feedback term given by 

poo 

F(r,x(f-r)) =7 / Gg(r)a;(f-r)dr (7) 

Jo 

with Gg('r) given by Eq. ([5]). For simplicity, we fix a = 1. 

A. Theoretical Estimation of the Response Amplitude 

We can determine the solution of ([ 6 ]) in the long time limit for D S> a; by writing 
x(t) = X{t) + fit) where X and V’ are the slow and the fast variables, respectively, and 
W = (l/ 2 vr)/o^” ^jJ dr = 0. The evolution equations for X and '0 are 

X + dX+ x + p{x^ + 

+F{t, X{t — t)) = f cosut, ( 8 a) 

'ip + d'lp + ujqiP + d^X'^ip + 3/3X — pip"^)) 

-f/5 (- 0 ^ — (' 0 ^)) + -^(u "Piflt — = g cos fit. ( 8 b) 
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As V’ is a rapidly varying function of time, it is reasonable to approximate Eq. ([ 8 ]d) as 

‘ijj + F{t,' ll!{Qt — Qt)) = g cos Qt. (9) 

For f —)■ oo, we assume the solution of Eq. ([9]) as Ah cos(f2f+ 0). The second term in Eq. ([9]) 
is then worked out as (for details see Appendix A) 

F{t, iliiytt — fir)) = Ah cos(fh: + 0 + 0), (10a) 


where 


p pn 

"-fj-T' 


(lOb) 


Substituting 0 = Ah cos{Qt + 0) and F given by Eq. fllOap in Eq. ([9]), we obtain 


Ah = —, /i = (fl^ — yfl ^cos0)^ + (7fl ^sin0)' 

/i L 


1/2 


= tan 


-1 


yfl ^sin0 
— 7 fl“P cos 6 


(lla) 

( ll b) 


Further, (0) = 0, (0^) = (l/27r) jQ^0^dr = Ah/ 2 and (0^) = 0. Then the equation of 
motion of the slow variable X given by Eq. ([ 8 ^) becomes 


X + dX + CiX + /3X^ + F{T,X{t - t)) = f cos ut, (12) 


where Ci = oOq + 30 ( 7 ^/( 2 /r^). 

Slow oscillations occur about the equilibrium points X* of Eq. flT^ with / = 0. To 
determine X*, we write X{t—r) = X{t) = X*. In this case F in Eq. (IT^ is simply 'yX* since 
/o°° = r(p). The equilibrium points are Xq = 0, and X^ = +>/—(Ci + 7)/0. 

For convenience, we introduce the change of variable Y = X — X* and obtain 

Y + dY + + 3(3X*Y'^ + (3Y^ + F{T,Y{t - t)) = f cos ut, (13) 


where a;/ = Ci + 3/3X*‘^. u,- is termed as the resonant frequency of the slow motion. For 
I/I -C 1 and l + l +; 1 we can drop the nonlinear terms in Eq. flT^ and write its solution in 
the limit f —)■ oo as + = Al cos{ut + <h). For this solution, referring to the Eqs. (lA.ip and 
flA.3p . the distributed delays feedback term F is written as 


F 


Re 


-7Al0(‘"*+'*’)' 
(l + ici;)P 


(14) 


5 

















In Eq. flA.Sp of Appendix A, the term 1/(1 + is approximated as i~PQ~P assuming 
that fl S> 1. We cannot use this kind of approximation in Eq. flT^ because cj can be < 1. 
However, defining a + i 6 = 1/(1 + we rewrite Eq. flT^ as 

F = 7 Al [acos(a;f + <!)) — 6 sin(a;f + $)]. (15) 

Following the procedure used to determine Ah, we obtain 

Al = //\/5, S = (cj^ - + (dcj + . (16) 

Then, we define the response amplitude Q = A^/f = l/\/S. 

B. Resonance Analysis 

We analyse the occurrence of a resonance in both the double-well and the single-well 
potentials of the system. The potential of the system in absence of damping, feedback term 
and the external periodic force is of a double-well form for cJq < 0 , /? > 0 and of a single-well 
form for cJq, /3 > 0 . 

We numerically compute Q, to verify the theoretical expression Q. For this purpose, we 
integrate Eq. ([ 6 ]) using the Euler method with step size 0.01. Leaving the solution corre¬ 
sponding to first 10 ^ drive cycles as a transient, we compute the sine and cosine components 


of Q, denoted as Qs and Qc, respectively. 

using the formula 


2 

pnT 


Qs 

nl 

/ x{t) sinuitdt, 

Jo 

(17a) 

2 

pnT 


Qc = ^ 
nl 

/ x{t) cosutdt, 

Jo 

(17b) 

where n = 10^ and T = ‘2'n ju. Then, Q = 

VQI + QVf- 



For the case of the double-well potential system, we fix the values of the parameters as 
d = 0.5, Uq = —1, /3 = 1, / = 0.1, w = 1, = 10 and p = 0.5. In Fig. [1^ both theoretically 

and numerically calculated Q are plotted as a function of the control parameter g for three 
fixed values of 7 . The theoretical Q closely matches with the numerically computed Q. 
For a wide range of values of g, the response amplitude for 7 = 0.1 is higher than that of 
< 5(7 = 0) , while Q (7 = —0.1) < (^(7 = 0). In Fig. [T^ Q becomes maximum at two values 
of g. In order to capture the mechanism of the observed vibrational resonance in Fig. [T]d, 
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FIG. 1: (a) Q versus g for different values of 7 for the system Q with a double-well potential. 
The continuous and dashed lines are theoretical and numerical results, respectively, (b) as a 
function of g. The horizontal dashed line corresponds to — 7 a. The two solid circles mark the 
values of g at which resonance occurs. The values of the parameters of the system ([^ are d = 0.5, 
Wq = —1, /3 = 1, / = 0.1, cj = 1, n = 10, and p = 0.5. 


the variation of = Ci + 3/3X*^ with g is shown for 7 = 0.1. For g < g^, where 


9c = 


2/i2 T 


, 7<|wol) 


(18) 


there are three equilibrium points = 0 and X^. Slow oscillations take place about X^. 
For g > gc there is only one equilibrium point Xq and a slow motion occurs about it. In the 
calculation of u/, we used X* = X^ = ±\/—(iFi + j)//3 for g < gc and X* = Xq = 0 for 


9 > 9c- 
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In the expression for Al given in Eq. flT 6 |) when g is varied, the quantity alone varies 
while all other terms remain the same. We notice that Q or Al becomes a maximum when 
the quantity S becomes a minimum. This happens whenever matches with — 7 a. In 
Fig. [T]d for 7 = 0.1 as 5 ^ increases from 0 the value of cj^ decreases from 2|a;g| — 87 = 1.7 
and becomes the minimum value —7 at g = gc- When g is further increased, increases 
monotonically. At two values of g, denoted as and g^^, — 7 a (marked by the 

horizontal dashed line in Fig. [T]d). At these values of g, the response amplitude attains the 
maximum value l/(duj + 76 ). For 7 = 0.1, the theoretical values of = 51 and g^^^ = 113, 
while the numerically predicted values are 47 and 112, respectively. Q becomes a minimum 
at g = gc at which is minimum and the number of equilibrium points changes from three 
to one. Essentially, at this value of g(= gc) the effective potential of the slow variable X 
changes from a double-well form to a single-well. 

From the theoretical expression of Q, it is possible to determine the values of g at which 
resonance occurs. From Eq. (1T6|1 we infer that Q becomes maximum when dS/dg = 0. This 
condition leads to the following results. 

(i) Resonance occurs at 


gd) = 

^VR 


r ,,2 


gd) = 


^ (2|i^ol + 7(0 - 3)) 

-|l/2 


11/2 




^ (lu,Sl + - 7a) 


> Qc 


provided 


UJ 


7ci = 


< 7 < 7c2 = 


2|a;ol — 


(19a) 
(19b) 

• ( 20 ) 

a — 1 3 — a 

This condition assures that Qc > 0, qd\ qd) > Q and qd^ < qc < 0 ^^^. The hrst resonance 
takes place at a value of g{= gflO) < g^ while the second resonance is at an another value 
fi'vR ^ fi'c- 9 = 9c the response amplitude attains a local minimum. The resonances 
and gd^ given by Eqs. (IT^ are due to the matching of with — 7 a. 

(ii) For 7 > 7 c 2 there is only one resonance at 5 ^ = gd^ = gc- Here, the maximum value 
ot Q at g = gc is due to the local minimization of the function S in Eq. 0161) for which 

-|- 7 a 7 ^ 0 . 

(iii) For 7 < 7 ci only one resonance is possible and it takes place at gd^. Fig. [2] presents the 
influence of the parameters 7 and g on Q for four hxed values of the parameter p appearing 
in the gamma distributed delays ( Eqs. ([7]) and ([5]). For p = 0.5, the conditions fl20|) are 
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(a) p=0.5 



(b) p=1.5 



(c) p=2.5 (d) p=3.5 




FIG. 2: Variation of Q as a function of 7 and g for four fixed values of p. The values of the other 
parameters are as in Fig. [TJ 

satisfied for 7 G [—4.49843,0.42857]. That is, two resonances can occur when g is varied 
from a small value for values of 7 in the above interval. Fig. [2^ conhrms this. For p = 3.5 
two resonance peaks can take place for 7 G [—0.78,0.31]. In Fig. |2l for all values of p the 
values of Q at resonance decrease with a decrease in the value of 7 from the value 1. As 7 
increases from 7 ci, both g^^^ and move away from gc. 

In Fig. the variation of 7 ci and ^2 with p is plotted. In the stripped regions two 
resonances occur. For the values of 7 and p above the curve 7 c 2 , there is only one resonance 
8-^9 = dyl > Qc given by Eq. flT^i while for below the curve 7 ci a single resonance occurs 
Sit g = gc- The Figs. [2]and|3^ clearly demonstrate the strong influence of the parameters of 
the delay feedback term. In Fig. |3 ]d theoretical and numerically computed g^-^ are shown for 
a range of values of 7 with p = 0.5. The theoretical g^-^^ matches closely with the numerical 
g^-^ for I7I < 1. For I7I > 1 the deviation between the theoretical g^-^ and the numerical g^-^ 
increases with an increase in | 7 |. 


9 






2 




7 


FIG. 3: (a) 7 ci and 7 c 2 versus p for the double-well potential Duffing oscillator. Two resonances 
take place in the stripped region while only one resonance occurs in the remaining regions, (b) 
versus 7 for p = 0.5. The continuous curve and the solid circles are theoretically and numerically 
computed values of respectively. 

Figure 0] depicts the change in the slow motion X{t) as a function of ( 7 , computed by 
solving Eq. (IT^ for the values of the parameters used in Fig. [T]with 7 = 0.1. For p = 0.5 
and 7 = 0.1, the values of and are 47, 78, and 112, respectively. For g < g^ 

there are two orbits with one centered at and another centered at X^ (not shown in 
Fig. 0]). The Q of both orbits are the same. In Fig. 01 the orbits numbered as 1, 2 and 3 are 
for g = 1, 25 and 47(= These are the orbits with X^ as the center. The orbits with 
X* as the center are not shown in Fig. 01 We can clearly see that Q of the orbit-3 is the 
maximum. Further, X^ (the center of the orbit) moves towards origin with increase in the 
value of g. At g = g^ = 78, the center of the orbit (orbit-4) becomes Xq =0. As g increases 
further, the center of the orbit remains the same, however, Q increases, reaches a maximum 
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FIG. 4: Phase portrait of the slow motion X{t) of the system ([ 6 ]) for different values of g with 
d = 0.5, ujq = —1, (3 = 1, f = 0.1, w = 1, n = 10, p = 0.5, and 7 = 0.1. The values of g for the 
orbits 1 — 6 are: g = 1(1), 25(2), 47(3), 78(4), 112(5), 185(6), respectively. 



FIG. 5: Variation of Q with oj in the absence of time-delayed feedback and high-frequency force 
for the single-well Duffing oscillator with d = 0.5, Wq = 1, /3 = 1 and / = 0.1. 


at g = = 112 (orbit-5) and then decreases (orbit -6 with g = 185). 

Next, consider the system ([6]) with a single-well potential (cUq, /5 > 0). In absence of delay 
and an external periodic force, the system flT^ has only one real equilibrium point Xq = 0 . 
The theoretical expression for g^-^ is given by 


9vk 


■2/r2 





( 21 ) 


At most one resonance is possible for —Wq — 7 a > 0 and no resonance for —cUq — 7 a < 0. 
Fig. Oshows Q versus a; for 7 = 0, 5 : = 0 with d = 0.5, Wq = 1, /? = 1 and / = 0.1. This figure 
shows the typical nonlinear resonance due to the applied single frequency periodic force. For 
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FIG. 6 : (a) Resonance region (stripped region) in the 7 versus p parameter space for the single-well 
Duffing oscillator with d = 0.5, Wq = 1, /3 = 1, / = 0.1, D = 10 and uj = 2. (b) Q versus g for four 
fixed values of 7 with p = 0.5. 

a range of values of w, the response amplitude Q is much lower than its value at resonance. 
For example, at ca = 2 the value oi Q = 0.31623 is lower than the value Q = 2.06527 
corresponding to a; = (Umax = 0.94. The value of Q at a; = 2 can be enhanced and resonance 
can also be realized by including the time-delayed feedback and the high-frequency force. 

From Eq. fl2T|) the condition for resonance for cUq = 1 and a; = 2 is 3 — 7 a > 0. If a > 0 
(a < 0), then a resonance occurs only for 7 < 7 ci = 3/a (7 > 7 c 2 = ~3/|a|). Fig. | 6 ^ depicts 
the plot of 7 ci and 7 c 2 . In the stripped region a resonance can takes place. For p = 0.5, the 
value of a is 0.56886 and hence a resonance can occur, when g is varied for 7 < 7 ci = 5.27371. 
In Fig. [6)d for 7 = 0 , and ±0.5 a resonance occurs. For 7 = 6 > 7 ci when g is increased 
from zero, the response amplitude monotonically decreases with an increase in g as shown 
in Fig. [6)d and there is no resonance. In the single-well case also a resonance occurs when 
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ul = — 7 a (Eq. ffT 6 |) ). where = Ci = Wg + j (2/i^). Unlike the double-well system 

the center of the slow motion (X(t)) and the actual motion {x{t)) always takes place about 
the origin. 


III. INTEGRATIVE TIME-DELAYED DUFFING OSCILLATOR 

In this section we present the results of the investigation of the vibrational resonance 
in the Duffing oscillator system with an integrative time-delayed feedback term given by 
Eq. (H]). 


A. Theoretical Response Amplitude 


We follow the theoretical procedure applied for the system with a gamma distributed 
time-delayed feedback with x(t) = X{t) + fit). We obtain the following results: 


Y{t) 

5 

A* 


Al cos(a;t -|- <h), 

— Q = — 

x/:s’ ^ x/:s’ 

2(2 7 . 

oj — \ u -sm ua 

V oja 


Ui + Ui = + 

0 , ±\/—(Ui -|- 7 )/13 , 


+ 

2 /i 2 ’ 


'1 

doo - (1 — coscaa) 

ooa 


OJ. 


fla 


^ -|- -Y— sin fla 


+ 


dfl — 7^(1 — cos Da) 
\la 


( 22 a) 

( 22 b) 

( 22 c) 

( 22 d) 

( 22 e) 

( 22 f) 


Notice the difference between the S^s given by Eqs. ffT 6 |) and fl22bi. Since the solution Y{t) 
is assumed to be periodic with period juj., we choose 0 < a < 2tt/ uj. 


B. Double-Well Potential System 


For cUq < 0, > 0 and 7 < |a;g| the values of g at which the number of equilibrium points 

changes from three to one and the resonance occurs are given by 

1/2 


gc = 


2 /i^ 


(l^ol -7) 


(23) 


13 














and 


— 

^VR 


^( 2 ) = 
^VR 


V 

3/3 

2 /i2 


( 2 |a;n| — 37 — o;^ + — sincua ) 
V ooa / 


11/2 


< 9c, 


2 1 2 'y ' 

Wn + cj-sincna 

ooa 


1/2 


> 9c- 


For two resonances to occnr the condition on 7 is jci < 7 < 7c2 where 

,2 ol, ,21 , ,2 


a; 


7ci = 


2|a;Q| — u" 


’ 1 1 • ’ 7c2 q 1 • 

1 -sin ooa 3 -sin oia 


( 24 ) 

( 26 ) 


(26) 


When 7 > 7 c 2 , only one resonance is possible with the corresponding given by For 
7 < 7 ci there will be only one resonance at 5 ^ = 

Fignre[7K depicts the threshold cnrves 7 ci and 7 c 2 for Uq = —1, (3 = 1, d = 0.5, / = 0.1, 
00 = 1 and 11 = 10. In the stripped region, two resonances occnr while only one resonance 
occnrs in the remaining region. In Fig. [Tb, both theoretically predicted and nnmerically 
compnted and are plotted as a fnnction of 7 for a = 1. Here again, the theoretical 
resnlt closely matches with the nnmerical simnlation. 

For a = 1, we find jd = —6.308 and 7 c 2 = 0.46394. Two resonances are possible for 
—6.308 < 7 < 0.46394 and only one resonance can occnr ontside this interval of 7 , when g 
is varied. When 7 = 0.2, is eqnal to — ('y sino;a)/ua at two valnes of g, namely at 
g = 43.5 and g = 113.55. At these two valnes of g, Q becomes maximnm. 00 “^ matches with 
— ( 7 sina;a)/a;a at only one valne of g for 7 = 0.5 and hence only one resonance. Here, 
the resonance occnrs at = 103.5 > gc = 59. 

The inflnence of the parameters g and 7 on Q is depicted in Fig. [H]for fonr fixed valnes 
of time-delay a. In this fignre, we clearly notice that for 7 valnes below a critical valne 
( 7 ci) only one resonance takes place. Comparing the Fig. [2] of the gamma distribntive time- 
delayed feedback case and Fig. [ 8 ] of the integrative time-delayed feedback case, the effect of 
a is fonnd to be similar to that of p. 


C. Single-well Potential System 


For the single-well potential case {oJq ,(3 > 0 ), the slow motion takes place abont the 
eqnilibrinm point X* = 0 and the theoretically predicted g^^^ is 
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2 /i 2 


2 2 7. 

—oon + u -sm ooa 

ooa 


1/2 


(27) 
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FIG. 7: (a) Dependence of 7 ci and jc 2 with the time-delay a for the system ([6j) with an integrative 
time-delayed feedback term. Here d = 0.5, Wq = —1, /3 = 1, / = 0.1, iv = 1 and D = 10. Two 
resonances occur in the stripped region. Below and above the stripped region only one resonance is 
possible, (b) Theoretically predicted (continuous curve) and numerically computed (solid circles) 
values of and for a = 1 as a function of the parameter 7. 

From (|?7|l . we find the condition for resonance as 

— -^sina;Q!> 0 . (28) 

ua 

We fix d = 0.5, cUq = 1, /? = 1, / = 0.1, uj = 2 and = 10. In this case, when g is varied Q 
can exhibit a resonance only if 

6 cr 

7 < 7 ci = . ^ , a e [0, 7 r/ 2 ] (29a) 

sin 2a 
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(a) a= 1 (b) a= 2 



(c) a= 3 



(d) a=4.5 

3r 



FIG. 8 : Q versus g and 7 for four fixed values of a for the double-well system ([ 6 ]) with an integrative 
time-delayed feedback term. 


or 


7 > 7c2 


—6a 

sin2Q! 


a e [7r/2, vr]. 


(29b) 


where the range of a is restricted to [0, 27r/a;] = [0, tt]. The threshold curves 7 ci and 7 c 2 are 
plotted in Fig. [9ti. In the stripped region, only one resonance can occur. In the remaining 
set of values of 7 and a, the integrative time-delay suppresses the existing resonance for 
7 = 0. When a = 0.2 < tt, the theoretically predicted 7ci = 3.08, while the numerically 
computed 7 ci = 3.04. In Fig. [9 )d both theoretically and numerically calculated Q are shown 
for four values of 7. For 7 = ±1 < 7ci, Q displays a single resonance. For 7 = 4 > 7^1, the 
response amplitude monotonically decreases with g and there is no resonance. 
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FIG. 9: (a) Single resonance region (stripped region) and no resonance region (blank region) for 
the single-well potential of the system ([ 6 ]) with an integrative time-delayed feedback term, yd and 
7 c 2 are given by Eqs. ([2^ . Here ooq = 1 and oj = 2. (b) Q versus g for four fixed values of 7 with 
a = 0 . 2 . The continuous and dashed lines are theoretically and numerically computed values of Q, 
respectively. 

IV. CONCLUSIONS 

In this paper, we have analyzed the influence of distributed time-delays and integrative 
time-delays on vibrational resonance for the Duffing oscillator, by using a theoretical ap¬ 
proach. Using the approximate analytical expression for the response amplitude, we have 
been able to determine the parameter regions (describing the time-delayed feedback) where 
either two resonances, one single resonance or no resonance occur. An enhanced response is 
realized for a range of values of the control parameters of the feedback. With and without 
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time-delayed feedback at most two resonances and at least one resonance occur in the case of 
the double-well Duffing oscillator. In the single-well case, at most one resonance is possible. 
Suppression of this single resonance by time-delayed feedback happens for a range of values 
of the parameters. Thus, the time-delayed feedback can be used to control the number of 
resonances and the value of the response amplitude at the resonance. 
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Appendix A: 


PROOF OF EQUATION IfTOah 

We present the details of the evaluation of the term — Dr)) given by Eq. ([7]) 

for -0 = Ah cos{flt + 0). Equation ([7]) with the "0 given above takes the form 


F = 


7^H 

r(p) 


^e cos(Dt + (j) — Qt) dr. 


Using cos 0 = Re (e'^) the above equation is rewritten as 


F = Re 


T^h , 

Lr(p) 


(nt+0) 


T 


e-’^ dr 


where 6 = 1 -I- iD. As the value of the integral in flA.2|) is r{p)/bP, we have 


F = Re 

= Re 


. (1 + iD)' 


i-PQ-P ( 1 - ^ 


-p 


For D > 1, (1 - ^) 


-p 


(g-i/O) P ^ gip/0_ 

F = Re 

= 'yAii^l~P cos{ilt + (j) + 9), 


where 0 = --—. 

D 2 


(A.l) 

(A.2) 

(A.3) 

(A.4) 
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